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Abstract 

Consider three qubits A, B, and C which may be entangled with each other. 
We show that there is a trade-off between A's entanglement with B and its 
entanglement with C. This relation is expressed in terms of a measure of 
entanglement called the "tangle," which is related to the entanglement of 
formation. Specifically, we show that the tangle between A and B, plus the 
tangle between A and C, cannot be greater than the tangle between A and 
the pair BC. This inequality is as strong as it could be, in the sense that for 
any values of the tangles satisfying the corresponding equality, one can find 
a quantum state consistent with those values. Further exploration of this 
result leads to a definition of the "three-way tangle" of the system, which is 
invariant under permutations of the qubits. 



PACS numbers: 03.65.Bz, 89.70. +c 



1 



Quantum entanglement has rightly been the subject of much study in 
recent years as a potential resource for communication and information pro- 
cessing. As with other resources such as free energy and information, one 
would like to have a quantitative theory of entanglement giving specific rules 
about how it can and cannot be manipulated; indeed such a theory has begun 
to be developed. The first step in building the theory has been to quantify 
entanglement itself. In the last few years a number of entanglement mea- 
sures for bipartite states have been introduced and analyzed [1—7], the one 
most relevant to the present work being the "entanglement of formation" 
M, which is intended to quantify the amount of quantum communication 
required to create a given state. In the present paper we draw on some of the 
earlier work on entanglement of formation || |7| in order to explore another 
basic quantitative question: To what extent does entanglement between two 
objects limit their entanglements with other objects? 

Unlike classical correlations, quantum entanglement cannot be freely shared 
among many objects. For example, given a triple of spin-| particles A, B 
and C, if particle A is fully entangled with particle B, e.g., if they are in the 
singlet state ^(|TI) — lit)); then particle A cannot be simultaneously entan- 
gled with particle C. (If A were entangled with C, then the pair AB would 
also be entangled with C and would therefore have a mixed-state density 
matrix, whereas the singlet state is pure.) One expects that a less extreme 
form of this restriction should also hold: if A is partly entangled with B, then 
A can have only a limited entanglement with C. The first goal of this paper 
is to verify this intuition and express it quantitatively. We will see that the 
restriction on the sharing of entanglement takes a particularly elegant form 
in terms of a measure of entanglement called the "tangle," which is closely 
related to the entanglement of formation. Further analysis of this result will 
lead us naturally to a quantity that measures a three-way entanglement of 
the system and is invariant under all permutations of the particles |3. 

The present work is related to recent work on the characterization of 
multiparticle states in terms of invariants under local transformations [9- 
12]; indeed, both the tangle and our measure of three-way entanglement are 
invariants in this sense. Our work is also related to research exploring the 
connection between entanglement and cloning [13-17]. An example along 
these lines was studied by Brufi, who asked, in the case of a singlet pair 
AB, to what extent particle B's entanglement with particle A can be shared 
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symmetrically and isotropically with a third particle, for a purpose such as 



teleportation where isotropy is desired [1£|. Our investigation is similar in 
spirit to that of Brufi but has a different focus in that we are looking for 
a general law governing the splitting or sharing of entanglement; thus, for 
example, we make no assumptions about the symmetry of the system. Some 
of the results presented here have been mentioned in a recent paper by one 



of us ||19|| , but the proofs and most of the details and observations have not 
been previously published. In this paper we confine our attention to binary 
quantum objects (qubits) such as spin-| particles — we will use the generic 
basis labels |0) and |1) rather than ||) and | j.) — but the same questions could 
be raised for larger objects. 

We begin by defining the tangle. Let A and B be a pair of qubits, and 
let the density matrix of the pair be Pab, which may be pure or mixed. We 
define the "spin-flipped" density matrix to be 

Pab = K ® (yy)p\ B {a y <g> a y ), (1) 

where the asterisk denotes complex conjugation in the standard basis {|00), 

1 1) , 1 10), 1 11)} and a y expressed in the same basis is the matrix ^ ^ 0^ ) 

As both pab and Pab are positive operators, it follows that the product 
Pab Pab, though non- Hermit ian, also has only real and non- negative eigen- 
values. Let the square roots of these eigenvalues, in decreasing order, be Aj, 
A 2 , A3, and A 4 . Then the tangle of the density matrix pab is defined as 

tab = [max{Ai - A 2 - A 3 - A 4 , 0}] 2 . (2) 

For the special case in which the state of AB is pure, the matrix PabPab has 
only one non- zero eigenvalue, and one can show that tab — 4detp J 4, where 
Pa is the density matrix of qubit A, that is, the trace of pab over qubit B. 

It is by no means obvious from the definition that the tangle is a measure 
of entanglement for mixed states. This interpretation comes from previous 
work, in which a specific connection is established between the tangle and 
the entanglement of formation of a pair of qubits Jj]]. For the purpose of this 
paper it is sufficient to note that r = corresponds to an unentangled state, 
r = 1 corresponds to a completely entangled state, and the entanglement of 
formation is a monotonically increasing function of r.[] (The earlier work did 

1 The entanglement of formation is given by E = + \\J\ — r), where h is the binary 
entropy function h{x) = —a; log a; — (1 — x) log(l — x). 



3 



not define the tangle per se but rather the "concurrence," which is simply the 
square root of the tangle. It is with some hesitation that we introduce the 
new term "tangle" here instead of speaking of the square of the concurrence. 
However, for our present purpose the tangle does seem to be the more natural 
measure to use, and the discussion is significantly simplified by introducing 
this term.) At present, the tangle is defined only for a pair of qubits, not for 
higher- dimensional systems. 

We now turn to the first problem of this paper: given a pure state of three 
qubits A, B, and C, how is the tangle between A and B related to the tangle 
between A and C? For this special case — a pure state of three qubits — the 
formula for the tangle simplifies: each pair of qubits, being entangled with 
only one other qubit in a joint pure state, is described by a density matrix 
having at most two non-zero eigenvalues. It follows that the product PabPab 
also has only two non-zero eigenvalues. We can use this fact and Eq. (0) to 
write the following inequality for the tangle tab between A and B. 

tab — (Ai — A2) 2 = Ai 2 + A2 2 — 2A1A2 

= Tt(pabPab) - 2AiA 2 < Tr (pabPab)- (3) 

Here pab is the density matrix of the pair AB, obtained from the original 
pure state by tracing over qubit C. Eq. (|3|) and the analogous equation for 
tag allow us to bound the sum tab + Tac : 

tab + t ac < T^(pabPab) + ^(pacPac)- (4) 

The next paragraph is devoted to evaluating the right-hand side of this in- 
equality. 

Let us express the pure state |£) of the three-qubit system in the standard 
basis {\ijk)}, where each index takes the values and 1: 

10 = J2 a ijk\ijk). (5) 

ijk 

In terms of the coefficients a^, we can write Tt(pabPab) as 

Ty(PabPAb) = ^ a ijka* mn k e rnrn>£nn'a* m > n > p aiijip£iii£j>j, (6) 

where e i = — e w = 1 and e o = en = and the sum is over all the in- 
dices. We now replace the product e nn i€jij with the equivalent expression 
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dnj'Sn'j — S n jS n 'j/, and in the first of the two resulting terms (that is, the one 
associated with 5 n j>5 n 'j) we perform a similar substitution for e m m'U'i- These 
substitutions directly give us 

T^(PabPab) = 2detp A - Tr(p|) + Tr(p^), (7) 

where pA, Pb, and pc are the 2x2 density matrices of the individual qubits. 
Because each of these matrices has unit trace, we can rewrite Eq. (|7|) as 

Tr (pabPab) = 2(det p A + det p B - det p c ) . (8) 

By symmetry we must also have 

Tr {pacPac) = 2(det p A + det p c - det p B ). (9) 

Summing these last two equations, we finally get a simple expression for the 
right-hand side of Eq. (|) , namely, 

Tt(PabPab) + Tt(pacPac) = 4detpA- (10) 

Eqs. ([|) and (|T0|) give us our first main result: 

tab + t A c < 4 det p A . (11) 

We can interpret the right-hand side of Eq. (|TT|) as follows. Thinking of 
the pair BC as a single object, it makes sense to speak of the tangle between 
qubit A and the pair BC, because, even though the state space of BC is 
four-dimensional, only two of those dimensions are necessary to express the 
pure state |£) of ABC. (The two necessary dimensions are those spanned by 
the two eigenstates of psc that have non-zero eigenvalues. That there are 
only two such eigenvalues follows from the fact that A is only a qubit.) We 
may thus treat A and BC, at least for this purpose, as a pair of qubits, and 
so the tangle is well defined. In fact this tangle ta(bc) is simply 4 det pa as 
we have mentioned before. We can therefore rewrite our result as 

tab + t~ac < Ta(bc)- (12) 

Informally, Eq. flT2| ) can be expressed as follows. Qubit A has a certain 
amount of entanglement with the pair BC. This amount bounds A's en- 
tanglement with qubits B and C taken individually, and the part of the 
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entanglement that is devoted to qubit B (as measured by the tangle) is not 
available to qubit C. 

We will say more shortly about the case of three qubits in a pure state, but 
at this point it is worth mentioning a generalization to mixed states. If ABC 
is in a mixed state p, then ta/bc) is n °t defined, because all four dimensions 
of BC might be involved, but we can define a related quantity t™^ c ^ via the 
following prescription. Consider all possible pure-state decompositions of the 
state p, that is, all sets {(il>i,Pi)} such that p = Y^iPil^i^il- F° r eac h of 
these decompositions, one can compute the average value of ta(bc)' 

(ta(bq) = J2Pi T MBC)(4>i)- (13) 

i 

The minimum of this average over all decompositions of p is what we take to 
be Ta(bc)^P)- The following analogue of Eq. ( |P2|) then holds for mixed states: 

TAB + r AC <r^ C y (14) 

To prove this, consider the pure states \ifji) belonging to an optimal decompo- 
sition of p, that is, a decomposition that minimizes (ta(bc))- We can write 
our basic inequality, Eq. (|i~2l) , for each such pure state and then average 
both sides of the inequality over the whole decomposition. The right-hand 
side of the resulting inequality is Ta(bc)^P)^ which is what we want on the 
right-hand side. On the left-hand side we have two terms: (i) the average 
of the tangle between A and B over a set of mixed states whose average is 
Pab (*-e., Trc(p)), and (ii) the average of the tangle between A and C over 
a set of mixed states whose average is pac- It is a fact that the tangle is a 
convex function on the set of density matrices.0 That is, the average of the 
tangles is greater than or equal to the tangle of the average. In this case the 
tangles of the averages are tab = t(pab) and tac = t(pac)- The sum of 
these two tangles must thus be less than or equal to Ta{bc)^P)i which is what 
we wanted to prove. 

Returning to the case of pure states, one may wonder how tight the 
inequality (pT2|) is. Could one find, for example, a more stringent bound of 

2 It follows from Ref. that the concurrence is a convex non-negative function on the 
set of density matrices for two qubits. The tangle, being the square of the concurrence, is 
therefore also convex. 
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the same form, based on a different measure of entanglement? To address 
this question, consider the following pure state of ABC: 

|0) = a|100) + /3|010) + 7|001), (15) 

where the three positions in the kets refer to qubits A, B, and C in that 
order. For this state, one finds that tab = 4|a| 2 |/?| 2 , tac = 4|a| 2 |7| 2 , and 
t a{ bc) = 4|a| 2 (|/?| 2 + M 2 ). Thus the inequality flT2|) becomes in this case an 
equality: tab + Tac — t a(bc)- This example shows that for any values of the 
tangles satisfying this equality, there is a quantum state that is consistent 
with those values. 

Now let T(r) be a monotonically increasing function of r that we might 
propose as an alternative measure of entanglement. For simplicity let us 
assume that T(0) = and T(l) = 1. Because of the above example, T could 
satisfy the inequality Tab + Tac < ^a(bc) only if T(x) + T(y) < T(x + y) for 
all non-negative x and y such that x + y < 1. Suppose T has this property. 
Then could there exist some quantum state for which Tab + Tac — ^a(bc) 
but tab + t ac < t a(bc)7 That is, could T yield an equality for some state for 
which t gives only an inequality? The answer is no, because if tab + t ac < 
ta(bc), then T A b + T A c = T(t A b) + T(t A c) < T(t A b + r A c) < T(t A {bc)) = 
Ta(bc)- Moreover, the only way T can match r in those cases where r gives 
an equality is for T to be equal to r. In this sense, r is an optimal measure 
of entanglement with respect to the inequality given in Eq. (0). Note, 
however, that the above argument applies only to functions of r. There 
could in principle be other measures of entanglement that are not functions 
of t that could make an equal claim to optimality. 

The entanglement of formation is a function of r, but it is a concave 
function and therefore does not satisfy an inequality of the form of Eq. (|T2|). 
Consider, for example, the state ^|100) + ||010) + ||001>. One finds that 
the relevant entanglements of formation are Eab = 0.601, Eac — 0.601, 
and Ea(bc) = 1- Thus, contrary to what one might expect, the sum of the 
entanglements of formation between A and the separate qubits B and C is 
greater than 1 "ebit," despite the fact that one can reasonably regard 1 ebit 
as the entanglement capacity of a single qubit. This is not a paradox; it 
simply shows us that entanglement of formation does not exhibit this par- 
ticular kind of additivity. (This sense of "additivity" should not be confused 
with the additivity of entanglement when one combines pairs to make larger 
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systems J?]]. It is not known whether entanglement of formation satisfies the 
latter notion of additivity.) 

We have just seen that there are some states for which the inequality 
( |i2|) becomes an equality. Of course there are other states for which the 
inequality is strict. As we will see, it turns out to be very interesting to 
consider the difference between the two sides of Eq. This difference can 
be thought of as the amount of entanglement between A and BC that cannot 
be accounted for by the entanglements of A with B and C separately. In the 
following paragraphs we refer to this quantity as the "residual tangle." 

Let the system ABC be in a pure state and as before, let the compo- 
nents of |£) in the standard basis be a^: 

10 = J2 a m\ i j k )- ( 16 ) 

ijk 

According to Eqs. @ and (|l7j|) and the discussion following Eq. ([LTD , the 
residual tangle is equal to 

Ta(bc) - tab - tac = 2(A^ S A^ B + A^A^ 7 ), (17) 

where Xf B and X^ B are the square roots of the two eigenvalues of PabPab, 
and \f c and \ AC are defined similarly. We now derive an explicit expression 
for the residual tangle in terms of the coefficients a^. 

We focus first on the product A^A^ 8 . This product can almost be 
interpreted as the square root of the determinant of PabPab- But PabPab 
is an operator acting on a four- dimensional space, and two of its eigenvalues 
are zero; so its determinant is also zero. However, if we consider the action 
of PabPab only on its range, then A^ B A^ B will be the square root of the 
determinant of this restricted transformation. 

The range of PabPab is spanned by the two vectors \vo) = J2ij <kjo\ij^) 
and |t>i) = J2ij a iji\ij^)- (These vectors also span the range of Pab-) To 
examine the action of PabPab on this subspace, we consider its effect on 

vectors of the form x\vf) + y\v2) = \ % ) . This effect is given by 
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where R is a 2x2 matrix. The product Xf\ 2 is the square root of the 
determinant of R. One finds that 

(19) 

where the sum is over all repeated indices. (We have ordered the factors 
so as to suggest the expression pabPab from which R is derived.) Taking 
the determinant of R involves somewhat tedious but straightforward algebra, 
with the following result: 

B \% B = VdetR = | d 1 - 2d 2 + 4d 3 | , (20) 

where 

j 2 2,2 2,2 2,2 2 

"1 — a 000 a lll ' a 001 a 110 ' a 010 a 101 ' a 100 a 011i 

d 2 = Oooo&maoiiflioo + flooo a m&ioi a oio + Ooooflm a iio a ooi (21) 

+ a on a ioo a ioi°oio + «oii a ioo a no a ooi + aioi a oio a no a ooi; 
d% = Oooofliio a ioifloii + ttmaooiOoio^ioo- 

We can get a mental picture of this expression by imagining the eight coef- 
ficients aijk attached to the corners of a cube. Then each term appearing in 
di, d<ii or efa is a product of four of the coefficients such that the "center 
of mass" of the four is at the center of the cube. Such configurations fall into 
three classes: those in which the four coefficients lie on a body diagonal and 
each one is used twice (<ii); those in which they lie on a diagonal plane (d 2 ), 
and those in which they lie on the vertices of a tetrahedron (d 3 ). Within 
each category, all the possible configurations are given the same weight. 

This picture immediately yields an interesting fact: the quantity \i B \ 2 B 
is invariant under permutations of the qubits. (A permutation of qubits 
corresponds to a reflection or rotation of the cube, but each di is invariant 
under such actions.) This means in particular that we need not carry out 
a separate calculation to find A^A^ , since we know we will get the same 
result. We can now therefore write down an expression for the residual tangle: 

ta(bc) - tab - r A c = 4|di - 2d 2 + 4d 3 |. (22) 

Note that the residual tangle does not depend on which qubit one takes 
as the "focus" of the construction. In our calculations we have focused on 
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entanglements with qubit A, but if we had chosen qubit B instead, we would 
have found that 



t~b(ca) - T B c ~ r B A = 4|di - 2d 2 + 4d 3 |. 



(23) 



The residual tangle thus represents a collective property of the three qubits 
that is unchanged by permutations; it is really a kind of three-way tangle. If 
we call this quantity tabci we can summarize the main results of this paper 
in the following equation. 



In words, the tangle of A with BC is the sum of its tangle with B, its tangle 
with C, and the essential three-way tangle of the triple. As an example, 
consider the Greenberger-Horne-Zeilinger state ^(|000) + 1 1 1 1 ) ) [p0[| . For 
this state the tangle of each qubit with the rest of the system is 1, the three- 
way tangle is also 1, and all the pairwise tangles are zero (the qubits in each 
pair are classically correlated but not entangled). Thus Eq. (p4|) in this case 
becomes 1 = + + 1. 

Finally, we note that the expression for tabc in terms of di, d 2 , and d% 
[Eq. Q22"D] can be rewritten, after a little more algebra, in a more standard 



where the sum is over all the indices. This form does not immediately reveal 
the invariance of tabc under permutations of the qubits, but the invariance 
is there nonetheless. 

It would be very interesting to know which of the results of this paper 
generalize to larger objects or to larger collections of objects. At this point 
it is not clear how one might begin to generalize this approach to qutrits or 
higher dimensional objects, because the spin-flip operation seems peculiar to 
qubits. On the other hand, it appears very likely that at least some of these 
results can be extended to larger collections of qubits. The one solid piece of 
evidence we can offer is the existence of a generalization of the state \<p) of 
Eq. (|i~5"D to n qubits: 



|0) = ai|100 . . .0) + a 2 |010 . . . 0) + a 3 |001 . . .0) + • • • + a„|000...1). (26) 



t~a(bc) — t ab + tac + tabc- 



(24) 



form: 



1~ABC - 2 }^ O-ijkO'i'j'mO'npk'O'n'p' 




(25) 
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One can show that for this state, the following equality holds. 

T12 + T13H h Ti n = Tl(23...n)) (27) 

where the qubits are now labeled by numbers rather than letters. We are 
willing to conjecture that the corresponding inequality, analogous to Eq. (|12|), 
is valid for all pure states of n qubits. 

This work was supported in part by the National Science Foundation's 
Research Experiences for Undergraduates program. WKW is grateful for the 
support and hospitality of the Isaac Newton Institute. We would also like to 
thank Tony Sudbery and Noah Linden for discussions that helped streamline 
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